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Abstract 

We consider a classical field configuration, corresponding to intersection of two 
domain walls in a supersymmetric model, where the field profile for two parallel walls at 
a finite separation is known explicitly. An approximation to the solution for intersecting 
walls is constructed for a small angle at the intersection. We find a finite effective length 
of the intersection region and also an energy, associated with the intersection. 



1 Introduction 



The problem of domain walls has a long history both in its general field aspect [|1|] and in a 
possible relevance to cosmology Recently the interest to properties of domain walls has 
been given a new boost in supersymmetric models, which naturally possess several degenerate 
vacua, thus allowing for a multitude fU of domain wall type configurations interpolating 
between those vacua. Furthermore it has been realised @, [|] that at least some of these walls 
have rather distinct properties under supersymmetry and in fact their field profiles satisfy 
first order differential equations in analogy with the Bogomol'nyi - Prasad - Sommerfeld 
(BPS) equations ||. For this reason such configurations are called "BPS-saturated walls", 
or "BPS walls". It was also noticed J?j that, quite typically, there may exist a continuous set 
of the BPS walls, all degenerate in energy, interpolating between the same pair of vacua. A 
further analysis M of such set in a specific model revealed that the solutions in the set can 
be interpreted as two elementary BPS walls parallel to each other at a finite separation, the 
distance between the walls being the continuous parameter labeling the solutions. Since all 
the configurations in the set have the same energy, equal to the sum of the energies of the 
elementary walls, one encounters here a remarkable situation where there is no 'potential' 
interaction between the elementary walls. (This property is protected by supersymmetry 
and holds in all orders of perturbation theory.) 

The existence of multiple vacua in supersymmetric models naturally invites a consider- 
ation of more complicated field configurations, than just two vacua separated by a domain 
wall, namely those with co-existing multiple domains of different degenerate vacua j3]. This 
leads to the problem of intersecting domain walls. It should be noted, that by far not every 
conceivable configuration of domains is stable, e.g. any intersection of domain walls in a 
one-field theory is unstable, and the conditions for stability of the intersections of the walls 
in multiple-field theories [Q] generally allow only a well defined set of intersection angles, de- 
pending on relation between the energy densities of the intersecting walls. As is discussed 
further in this paper, for the elementary BPS walls, considered in |§, i.e. non- interacting in 
the parallel configuration (zero intersection angle (3), at least a finite range of values for (3, 
including /3 — 0, is allowed by the stability conditions for a quadruple intersection, shown 
in Fig.|l|. When viewed as an intersection of world surfaces of the walls in the space-time, 
rather than as a static spatial configuration, the intersection describes scattering of moving 
walls, and in the case of a collision of spatially parallel walls (or kinks in a (1+1) dimensional 
model) the angle /3 translates into (3 = v/c with v being the relative velocity of the walls, 
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for small (3, i.e. for non-relativistic collisions. The purpose of the present paper is a more 
detailed study of the field profile for the intersection configuration of such type. We consider 
here the same super symmetric model as in Ref.|| and use the explicit solution for parallel 
walls, found there, to construct an approximation to the classical field profile for intersecting 
walls, which is valid to the first order in f3 for small j3. This approximation gives the energy 
of the static field configuration up to (3 2 inclusive, while for the case of collision of the walls 
the corresponding result is obviously the action for the collision "trajectory" . 




Figure 1: Quadruple vertex for four different vacuum domains 



Our main findings are as follows. There exists a well defined gap Ax between the vertices 
of pairwise 'meet' of the walls, so that a more detailed picture of the intersection looks as 
shown in actual detail in Fig.|2| further in the paper. In the collision kinematics this gap is in 
time and corresponds to a finite time delay in the scattering process. This classical quantity 
clearly can also be translated into a phase shift in a quantum mechanical description of the 
scattering. The gap depends on the intersection angle, and is given by 

= (1) 

p m 

where m is a mass parameter for the masses of quanta in the model, and l§ is a dimensionless 
constant, depending on the ratio of the coupling constants in the model. We also find that 
there is a finite angle-dependent energy (in the static configuration) associated with the 
intersection: 

e = Pl*-, (2) 
m 

where \i is the energy density of each of the elementary walls. Naturally, in a (3+1) dimen- 
sional theory e is in fact the energy per unit length of the intersection, while in a (2+1) 
dimensional case it is an energy localized at the intersection. In a (1+1) dimensional model 
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only the collision kinematics (for two elementary kinks) is possible, thus e has the meaning 
of a finite action (phase shift) associated with the collision. We believe that the equations 
(0) and (0) are quite general at small /3, and the only model dependence is encoded in the 
dimensionless quantity £ . 

The rest of the paper is organized as follows. In Section 2 we present the supersymmetric 
model under consideration and describe the solutions for the BPS walls. In Section 3 we 
consider intersecting domain walls within an Ansatz allowing us to construct the field profiles 
in the limit of small intersection angle, and we find the characteristics of the intersection. In 
Section 4 we discuss applicability of the findings of this paper when some of the restrictions, 
assumed here, are relaxed. 

2 Domain walls in a SUSY model 

The specific SUSY model under consideration is that of two chiral superfields $ and X with 
the superpotential 

W{$, X) = — $--A$ 3 -a$X 2 . (3) 
A 3 

Here m is a mass parameter and A and a are coupling constants. The phases of the fields 
and of the W are assumed to be adjusted in such a way that all the parameters are real and 
positive, and the lowest components of the superfields $ and X are denoted correspondingly 
as <p and x throughout this paper. The model has the Z 2 x Z 2 symmetry under independent 
flip of the sign of either of the fields: $ — > — $, X — > —X. The vacuum states in this model 
are found as stationary points of the superpotential function (|3]): dW(<p, x)/9<p — and 
dW(cj), X)/dx = 0, and are located at = ±m/A, x = (labeled here as the vacua 1 and 
2), and at = 0, x = iwi/ VAa (the vacua 3 and 4). The locations and the labeling of the 
vacuum states are shown in FigfJ. Throughout this paper the fermionic superpartners of the 
bosons are irrelevant and also only real components of the fields appear in the considered 
configurations. Therefore for what follows it is appropriate to write the expression for the 
part of the Lagrangian describing the real parts of the scalar fields: 

L = {d<pf + {dx f - - A 2 - a x 2 ") 2 - 4a 2 2 X 2 . (4) 

As discussed in detail in Ref. f8|, in this model exists a continuous set of BPS walls 
interpolating between the vacua 1 and 2, all having the energy density /X12 = § ^jt- These 



3 



4> 

O 2 



Figure 2: The diagram showing locations of the four vacua and their labeling in the SUSY 
model with the superpotential of eq.(|3]). 



configurations with positive (negative) x can be interpreted as parallel elementary BPS walls, 
i.e. connecting the vacua 1 and 3 (1 and 4) and the vacua 3 and 2 (4 and 2) located at a 
finite distance from each other. The energy p, of each of the elementary walls is p, — /X12/2. 
The only non-BPS domain wall in this model is the one connecting the vacua 3 and 4, and 
its energy is given by J X/a The first-order equations for the BPS walls are conveniently 



written in terms of dimensionless field variables / and h, defined as 



m 



./ 



m 



X 



h . 



X " v VAa 

For a domain wall perpendicular to the z axis the BPS equations read as| 



dz 

dh 
dJ 



f-h* 



2 -fh 

p 



(5) 



Here the notation is used p = X/a and the mass parameter m is set to one. 

Although the equations (|5]) are solved |§ in quadratures for arbitrary p, it is only for 
p = 4 that the non-trivial solution can be written in terms of elementary functions. It is this 
explicitly solvable case of p = 4 that we consider for the most part in this paper for the sake 
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of presenting closed expressions wherever possible. The set of solutions to the equations (||) 
in this case reads as [[J 

/(*)= a }T~ l \ z i h\z)= 2 f - , (6) 
a + 2e z + ae zz a + 2e z + ae zz 

where a is a continuous parameter, < a < oo labeling the solutions in the set, and the 
overall translational freedom is fixed here by centering the configuration at z = in the 
sense that f(z = 0) = 0. The interpretation of these configurations becomes transparent, if 
one introduces the notation: 

cosh s = a -1 . (7) 

Then the expressions @ can be written as 

„, , 1 / , Z — S , Z + ,9/ x 1 / , z — s , z + s\ , , 

f(z) = - (tanh + tanh -^-J , h\z) = - (l - tanh tanh -^-J . (8) 

In this form it is clear that at large s (i.e. at small a) the functions / and h differ from their 
values in one of the vacua only near z = ±s, i.e. the configuration splits into the elementary 
walls, corresponding to the transition between the vacua 1 — > 3 ((/, h) = (—1, 0) — > (/, h) = 
(0, 1)) at z = — s and the transition between the vacua 3 —>■ 2 at z = s. (For definiteness we 
refer to the solution with positive h(z).) Remarkably, the function / is simply a sum of the 
profiles for the elementary walls: f(z) = f\${z + s) + f^{z — s) with 

i ftanh* - - • w • z 



Mz) = j l^tanh - - 1 j , f 31 (z) = - (^tanh - + . (9) 

while the function /i(z) is not a linear superposition of the profiles of h for elementary walls: 

h 2 13 (z) = \ (l + tanh |) , /&(*) = i (l - tanh , (10) 

but rather for large s the profile of h(z) is exponentially (in s) close to a linear superposition 
of the elementary wall profiles separated by the distance 2s. Thus at large s (small a) the 
solution ([EJ) describes two far separated elementary walls. 

The caveat of parameterizing the solution in terms of s for arbitrary positive a is that, 
according to eq. @ , at a — 1 the parameter s bifurcates into the complex plane and becomes 
purely imaginary, reaching s = ±iir/2 at a = 0. (The functions f(z) and h(z) obviously are 
still real at < a < 

1 It can be noted that the profile of the fields at a = 1 is given by the special solution to the BPS equations 
first found in Ref.Q] at any p, such that p > 2. 
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3 Intersecting domain walls 



In order to describe a static intersection of two elementary walls at a small angle /3, we 
make the natural Ansatz that the parameter a in the solution (^) is a "slow" function of the 
coordinate x, i.e. / — > f(a(x), z) and h — > h(a(x), z). We then substitute this Ansatz in the 
expression for the energy of the fields and find 



E 



dx 




with 



F 2 (a) 



1 



(a 2 - 1 



5 3 a 2 arctan y? 
2 + 2 v / ^ r T 



1 



(12) 



The interpretation of this expression for the energy becomes quite simple at small a, if one 
writes it in terms of the (x dependent) parameter s (cf. eq.(|7|)), assuming that a < 1: 

2" 



E 



dx 



1 / , 2 5 3s 

2 H 5— cosh 2 s 1 — 

sinh s V 2 2 smh s cosh s 




(13) 



At large s the weight function for (ds/dx) 2 in the latter expression rapidly reaches one, and 
for the trajectory s = flx/2, corresponding to two walls being at large distance and inclined 
towards each other at the relative sloped /3, the energy per unit length of x coincides in order 



(3 with that of two independent walls: 2 fi Jl + (f3/2) 2 . 

In order to find the relation between a and x at arbitrary separation between the walls 
within our Ansatz, one needs to solve the variational problem for the energy integral, given 
by eq . (|TT|) . The solution is quite straightforward: introduce a function a(x), such that 

da 



da 

dx 



-F(a) 



dx 



(14) 



(The minus sign here ensures that a is growing when a decreases.) Then the solution of the 
variational problem for a is a linear function of x: 



a 



px/2 . 



(15) 



The overall shift in x and a is chosen so that o = and x = when a = oo i.e. at the 

center of the intersection. The slope of the cr(x) is determined by noticing that at a — > the 

2 At small j3 we make no distinction between the angle and the slope. When discussing large f3, or 
higher order terms, we imply that /3 is the relative slope, so that the full opening angle between the walls is 
2 arctan(/3/2). 
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function F(a) behaves as F(a) = a~ l + 0(a), thus at small a the slope of a coincides with 
that of s(x), defined above as (3/2 at large s. 

At this point we address the question about the accuracy of our Ansatz at small (3, and 
consider the full second-order differential equations for the fields <ft and \ following from the 
Lagrangian in eq.(£|). Since at fixed a the profile given by eqs.(^) satisfies also the second- 
order equations in the z variable, the mismatch in the full two-dimensional equations is given 
by the second derivatives in x only: d 2 f(a(x), z)/dx 2 and d 2 h(a(x), z)/dx 2 . One can readily 
see however that these quantities are of order f3 2 . Indeed, e.g. for the term with f(a(x),z) 
one finds using eqs. flliD and fllSf) 

d 2 f(a(x),z) = (3 2 ( d 2 f(a,z) _ F\a) df(a,z) \ 
dx 2 4F 2 (a) V da 2 F(a) da ) ' 

Thus our Ansatz correctly approximates the actual solution in the first order in (3. Once it 
is established that the correction (fi^ to the solution within the Ansatz is 0((3 2 ) (with 0j 
generically denoting the fields <p an d x)i if is clear that the corrections to the energy can 
start only in the order f3 A . Indeed, the variation of the energy, minimized within the Ansatz, 
SE/S(pi is of order of the mismatch in the field equations, i.e. 0(f3 2 ). Thus the error in the 
found energy 

A remark concerning further details of the discussed solution is due in relation with 
the singularity in a at x = 0. Formally, the described solution is only specified so far at 
x > 0. At negative x one a priori can choose one of two options: symmetrically reflect the 
profile at positive x, or also flip the sign of the field h at negative x. The first option would 
correspond to the domains of the same vacuum (i.e. the vacuum 3 or the vacuum 4) at small 
z on both sides of the intersection, while the second option describes the change from the 
vacuum 3 to the vacuum 4 at the intersection. The first configuration however is unstable 
(the translational zero mode develops a nodal line), and only the second type configuration 
should be chosen. 

One can notice that the profile of the fields within our Ansatz depends in fact on the 
scaling variable a = (3 x/2. This obviously implies that a 'longitudinal' interval of Sx between 
some fixed characteristic values of the fields in the configuration scales as (3~ l . The behavior 
of the parameter a in the expression (|6]) for the field profile as a function of a is determined 
by the equation 

„ ; , da . 
F(a) — = -1 . (16) 

This equation can be readily solved numerically, and the resulting picture of the intersection 
is presented in Figures 3 and 4. 
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Figure 3: Contour plot of the energy density at the intersection of the domain walls as de- 
scribed within our Ansatz in terms of the scaling variable (5 x/2. Darker shading corresponds 
to larger energy density. 



It is seen from the figures that the four crests of the energy density profile do not intersect 
at one point, rather there is a finite gap between the pairwise intersections of the energy 
crests for the 'initial' and the 'final' walls. In order to quantify this gap, we find from eqs.(|l6l) 
and (0) the relation between a and s as 



a = / F(a) da (17) 

J 1/ cosh s 

and then find numerically that at large s this relation gives a = s + 5 with 5 = 1.7527. . .. 
Since at large s the distance between the elementary walls is naturally identified as 2s, one 
concludes that when the intersection of two walls is extrapolated from large distances it 
comes short of the actual center of the solution by the distance 25 / (3. The gap £o between 
the apparent pairwise collision vertices is then twice this distance. In this way we obtain the 
equation ([J) with £q = 4 5 = 7.0108 ... in the particular model with p = 4. 

In order to find the proper expression for the energy associated with the intersection, 
one can consider the problem of intersecting walls as a boundary value problem in a large 



Figure 4: Three dimensional plot of the energy density at the intersection. 



box with the size L in the x direction (so that the boundaries are at x = ±L/2), and 
at the boundaries the positions of the walls in the z direction are specified as s = ±1/2, 
where s is the parameter in the expression (J8|) for the profile. (Clearly the sign of h(z) at 
x = +L/2 should be chosen opposite to that at x — —L/2 in order to have an intersection 
configuration.) We also assume that the separation I between the walls at the boundaries 
satisfies I <C L in order to correspond to /3 <C 1. For free walls, with no energy associated 
with the intersection, the energy of such configuration would be determined by the total 
length of the walls in the box: 

E = 2/i VL 2 + P w 2/iL + fiL (jj . (18) 
For the solution within our Ansatz, however, the energy is given by 

E = 2/xL + fiL^j , (19) 

where j3/2 is the slope in the solution a(x) = (3x/2. At x = L/2 one can use the asymptotic 
relation a = s + 5, and find from the boundary condition for s the relation for (3: 

fL-L — 

2~L + ~L' 
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Upon substituting this relation in eq. ([19]) one finds a finite difference between the expressions 
(|I9|) and ([18]) for the energy: e = E — E Q , given by our result in eq.(0). 

If instead of the spatial coordinate x one uses the time coordinate t, the discussed config- 
uration describes a collision of two parallel walls (or two kinks in (1+1) dimensional theory). 
In this case the parameter of the intersection (5 is identified as the relative velocity v of the 
walls, and eq.(|l]) describes the time delay in the process of scattering of the walls, while 
the quantity e in eq.(||) gives the additional action with respect to the free motion of the 
walls, and thus is equal to the phase of the transmission amplitude in a quantum mechanical 
description of the scattering. 

4 Discussion 

The description of intersection of elementary domain walls is found here for small intersection 
angles and for a particular value of the parameter p, p = 4, in a specific SUSY model. 
At present we can only speculate how our results are modified beyond these restrictions. 
Although the same Ansatz, as used here, can be applied in the limit of small f3 at any value 
of p, one inevitably has to resort to numerical analysis because of lack of an explicit simple 
expressions for the profiles of the fields, except for the trivial case p — 1. In particular, the 
existence of the gap at the intersection appears to be related to the existence of the special 
bifurcation solution |7J for parallel walls. Indeed a non-zero difference S between a and s, 
arising from the integral in eq. (|TT|) , comes mostly from the integration over the values of a 
past the bifurcation point, i.e. from a = 1 to a = oo. The bifurcation solution on the other 
hand exists in the considered model only for p > 2, and at p = 2 the solution corresponds 
to x( z ) — 0, so that the bifurcation point is located at x = for such p. Our preliminary 
numerical study of the dependence of the size of the gap on p indicates that indeed the gap 
appears to arise starting with p = 2 and becomes larger with increasing p, whereas at p < 2 
we have found no apparent gap. This conjecture about such behavior is somewhat supported 
by the fact that at p — 1 the elementary walls are not interacting and simply "go through" 
at the intersection with no gap whatsoever. 

As to the dependence on the intersection parameter j3, our Ansatz becomes inapplicable 
when f3 is not small, and the behavior of the corresponding configurations is not known. Here 
we can only mention that the stability conditions for the wall intersections allow a static 
quadruple intersection with any j3 if p > 1. Indeed in this case the energy of the non-BPS 
wall (3 — > 4), /i34 = 2 y/p p is larger than the sum of the energies of the elementary walls, and 
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thus the quadruple intersection cannot split into triple ones, since the stability conditions for 
triple intersections can not be satisfied. On the contrary, at p < 1 such splitting is possible 
with the critical value of (3 defined as f3 c / (2 yjl + (/3 c /2) 2 ) = ^fp. 

BPS solitons are known to present effective low-energy degrees of freedom in some super- 
symmetric field theories. These effective ("dual") theories exploit the fact that elementary 
solitons do not interact with each other at zero energies, both classically and quantum 
mechanically. This is just the case in the model considered, where the energy of parallel 
elementary domain walls does not depend on distance between them, and this degeneracy 
is not lifted by quantum corrections. We demonstrated here by an explicit calculation that 
some highly nontrivial interaction between BPS domain walls arises at nonzero momenta. If 
this observation holds for general case, to obtain dynamical information from dual models 
of that class appears to be an extremely complicated task. 
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